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Abstract 

Spontaneous symmetry breaking is an essential feature of modern science. We demonstrate that 
it also plays an important role in the physics of complex plasmas. Complex plasmas can serve as a 
powerful tool for observing and studying discrete types of symmetry and disordering at the kinetic 
level that numerous many-body systems exhibit. 
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I. INTRODUCTION 



In this Letter we address dynamical processes in highly ordered complex plasmas associ- 
ated with spontaneous symmetry breaking. 

Spontaneous Symmetry Breaking (SSB) plays a crucial role in elementary particle physics 
but is also very common in classical physics It happens whenever the system goes from 
a state which has a certain symmetry, e.g. rotational symmetry, into an ordered state, 
which does not have this symmetry anymore. In general, this state not necessarily has 
to be the ground (vacuum) state and the transition to the new state may or may not be 
associated with a phase transition. For example, in the case of magnetization the spins 
point all in one direction (ordered state) whereas above the Curie temperature there is no 
preferred direction. Another example from a mechanical system without phase transition is 
a vertical stick which bends under a sufficiently high force from above to one side breaking 
the rotational symmetry of the system without the force. 

Different symmetries coexisting in the same phase, and symmetry transformations escort- 
ing phase transitions are widely spread in nature. For instance, the mechanisms of symmetry 
breaking are thought to be inherent in the molecular basis of life [2]. SSB is also an im- 
portant feature of elementary particle physics [3j. The Universe itself is believed to have 
experienced a cascade of symmetry-breaking phase transitions which broke the symmetry of 
the originally unified interaction giving rise to all known fundamental forces 4|t6]. 

Symmetry effects are crucial either in 3D and 2D systems. Chiral (mirror-isomeric) 
clusters pi , magic clusters of a new symmetry 'frozen-in' by a solid surface [8( , or dynamical 
symmetry breaking by the surface stress anisotropy of a two-phase monolayer on an elastic 
substrate [<J are examples of the importance of 2D or quasi-2D systems in many applications. 

Low pressure, low temperature plasmas are called complex plasmas if they contain mi- 
croparticles as an additional thermodynamically active component. In the size domain of 
1-10 /im (normally used in experiments with complex plasmas) these particles can be vi- 
sualized individually, providing hence an atomistic (kinetic) level of investigations 101. 111]. 
The interparticle spacing can be of the order of 0.1-1 mm and characteristic time-scales are 
of the order of 0.01-0.1 s. These unique characteristics allow to investigate the microscopic 
mechanism of SSB and phase transitions at the kinetic level. 

Common wisdom dictates that symmetry breaking is an inherent attribute of systems in 



2 



an active state. Hence these effects are naturally important in complex plasmas where the 
particle cloud-plasma feedback mechanisms underlying many dynamical processes are easy to 
vitalize. Also in complex plasmas where different kind of phase transitions exist, e.g. in the 



elelectrorheological plasmas 



121 ] . one can find examples for classical SSB. Another option, 



interesting in many applications, is the clustering of a new phase which is dissymmetric 
with regard to a background symmetry (as an example of fluid phase separation in binary 



complex plasmas see 



13|). 



It is important to mention that the microparticles, collectin g el ectrons and ions from 
the plasma background, become charged (most often negatively [l0(]) and hence should be 
confined by external electric fields. The configuration of the confining forces might deeply 
affect the geometry and actual structure of the microparticle cloud. In rf discharge complex 
plasmas the particles are self-trapped inside the plasma because of a favorable configura- 
tion of the electric fields jl4]. One of the interesting things is the possibility to levitate a 
monolayer of particles under gravity conditions. In this case the particle suspension has a 
flat practically two dimensional structure. This is, of course, a very attractive simplification 
('from a theoretical point of view'), significantly lowering the description difficulties. Below 
we concentrate mostly on 2D complex plasmas. 

Depending on the discharge conditions, the monolayer can have crystalline or liquid order. 
2D configurations of dust particles either in crystalline or liquid state were successfully used 
to study phase transitions, dynamics of waves and many transport phenomena in complex 



plasmas 15l-l21|. A symmetry disorderin g es corting a crystalline-liquid phase transition has 



been investigated experimentally in 15Hl7j. Dislocation nucleation (a 'shear instability') 
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231 ] . albeit the importance of SSB for this phenomenon has not 



has been reported in 
been explained. 

The results of these recent experimental observations can not be properly addressed with- 
out a deep understanding of this important issue. We would like to highlight this in the paper 
and report on the physics of spontaneous disordering of a 'cold' plasma crystal, simulated 
melting and crystallization process, including associated defect clusters nucleation, disso- 
ciation, and symmetry alternation. These options are realizable in experimental complex 
plasmas, and can be mimicked in simulations, as we demonstrate below. 
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II. SPONTANEOUS DISORDERING OF A 2D PLASMA CRYSTAL 



It is well known that two broken symmetries distinguish the crystalline state from the 
liquid: the broken translational order and the broken orientational order. In two dimensions 
for ordinary crystals it is also well known that even at low temperatures the translational or- 
der is broken by spontaneous disordering mediated by thermal fluctuations 24|. As a result, 
the fluctuation deflections (disordering) grow with distance and translational correlations 
decay (algebraically, see [251] ). 

2D plasma crystals also obey this common rule. The character of disordering may be 
deeply affected by the confinement forces, though. Usually such an 'in-plane' confinement is 
due to the bowl-shaped potential well self-maintained inside the discharge chamber, which 



that is U, 



conf 



\MVL 2 r 2 , where 



to first order is approximately parabolic (see, e.g. [26 
r is the distance, M is the particle mass, and Q is the confinement parameter [28|. (The 
'out-of-plane' confining forces, controlling the position of the entire lattice, are normally 
much stronger; below we consider the 'pure' 2D-case, assuming, hence, an absolutely stiff 
out-of-plane confinement.) 

The fluctuation spectra can be calculated in the following manner. The long-range phonon 
contribution to the free energy of a 2D system of particles interacting via the Yukawa 
potential and confined by a shallow isotropic parabolic well can be conveniently represented 
as Q: 

M 



AU 



cm 



'1,2 



~ 2 4- fir 2 



5c 2 = n 2 /k 2 , 



(i) 

(2) 



where Ct r ,i are the transverse (shear wave) and the longitudinal (compressional wave) sound 
speed, and Vk, -Dk are the Fourier components of the vorticity V = curl z u and the divergency 
D = div u of the particle displacements u = Uk exp(zkr), and k is the wave vector 
(k = |k|). The unperturbed crystal is supposed to be hexagonal. 

The relationship ([1]) provides (see, e.g. 25j) the probability of the fluctuation w ~ 
exp(— AU/T). Next, using it, we can calculate the averaged fluctuation spectral intensity 
per unit mass as 



< u k >= 



J T 



+ 



k 2 c 2 tr + n 2 1 k 2 c 2 + n 2 ' T 

where vt is the particle thermal velocity. 



T 

M' 



(3) 



4 



FIG. 1: (color online) Cluster of 721 particles, interacting via the Yukawa type forces, confined 



inside the parabolic well (for simulation details see 281]) . The cluster size is R/a ~ 40. Shown are 
(to the left) the triangulated particle positions (defect 7- and 5-fold cells are colored in red and 
green) and (to the right) the cluster density profile. On average, the cluster density systematically 
decreases toward the edge. The local density maxima (minima) are easily identified as the positions 
of 5-fold (7-fold) cells. Note that the first circular row of the defects happens to appear at 0.46-R 
from the center. This agrees with the theoretical estimate ~ (0.44 — 0.48)i? obtained by using ([8]). 
Note also that the left part and the right part of the simulated cluster image are perfectly chiral 
(mirror-symmetric) . 



It is known (see 28|, |29|) that a lattice layer of a finite size R is stably confined if roughly: 



c; 



= p ~ const = 3 — 4. (4) 



Vt 2 Ra 

Here the parameter p stands for an effective number of the nearest neighbors of any edge 
particle. Note that according to (j4]) formally Q — > at R — > oc. 

Without confinement (f2 = 0) the fluctuation spectrum ([3]) apparently diverges oc k~ 2 
at k — > 0, and, as a consequence, in agreement with j^, 25], the crystal ordering decays 



algebraically with the distance r, i.e. the density-density correlation behaves as 

< p(ri)p(r 2 ) > -p 2 oc r~ n cos(br). (5) 

Here r = ri — r2, b is the vector of the reciprocal lattice, and n = It is assumed that 

r ^> a is large compared to the interparticle separation a. 
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In the experiments fl ^ is always finite (though noticeably small, one or two orders of 
magnitude less than the frequency of the local 'caged' oscillations of the individual particles 
27]). From at non- vanishing Q it immediately follows that the fluctuations remain 
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finite even at k — > 0. This absence of a singularity alters the character of disordering from 
algebraic (J5J) to exponential at a scale depending on the confinement parameter: 



r ~ r r 



ci/fl < R. 



(6) 



It is essential that both asymptotes - algebraic and exponential - must be treated as 'near- 
field' (r <C r c ) and 'far-field' (r ~ r c ) approximations. Hence it would be logical to assume 
that the ordering decay alternates with distance from algebraic to exponential. This is 
indeed in qualitative agreement with observations 30h32|. 



Remarkably ([I])-© are formally similar to the equations describing director fluctuations 



in nematic crystals in the presence of a magnetic field [25|, |33| . The action of the magnetic 
field is known as suppressing the large-scale director fluctuations in liquid crystals. 

The length scale ~ r c seems to be of a fundamental importance. The particles, experienc- 
ing a horizontal confinement, are distributed non-uniformly. The steady-state displacements 



of the particles u in the plasma crystal from their ideal 
represent a growing function with distance u = r 3 /8r^ 

lr 2 



u 



> L, 



ocations in a uniform 2D crystal 
The lattice breaks up when 

(7) 



where L is the Lindemann parameter. Since L=0.16-0.18 (see, e.g. 20[), it follows that the 
first row of defects most probably appears at r ~ r c ySL ~ (1.1 — 1.2)r c . Making use of fl4j), 
(]6]) one can estimate the size of the domains (or equivalent correlation length) as: 



R 



a 

S P L-. 



(8) 



The correlation length (JS} does not depend explicitly on the temperature. In other words, 
for purely topological reasons the big crystal spontaneously splits, assembling an array of 
sub-domains, even at zero temperature. The estimated values of r cor agree well with those 
obtained in the simulation - see Fig. [U and in experiments. For instance, it has been 



observed in 
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3l| that the crystal orientational order had a power law decay at distances 



r/R < 1/4 in fairly good agreement with w 0.3 following from (jSJ). 

The 'one-plus' correlation length ([8]), unavoidably introducing a network of sub-domains 
to a lattice layer, is of crucial importance, e.g., for observations of the so called hexatic state 



in the p 
studies 



.asma crystals that is still an outstanding and controversial issue in complex plasma 
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III. DEFECT CLUSTERS IN PLASMA CRYSTALS 



One of the possible scenarios for melting (recrystallization) in a 2D complex plasma 



is a precipitous increase (decrease) in the density of the dis 
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ocations and the dislocation 



34(. To realize this scenario 



aggregates (such as defect clusters, grain boundaries etc.) 
in simulations, it is desirable to avoid any aforementioned complications associated with the 
lattice layer sectioning 'from the very beginning'. A promising tool in that sense, allowing 
to create a defect-free initial lattice layer, is a hexagonal confinement cell proposed in 16]. 

We performed a series of simulations that revealed several peculiarities in symmetry that 
are worth to mention. 




<K>~\ eV" 1 



FIG. 2: Dislocation density vs. inverse mean particle kinetic energy. The dashed line is the 
least squares fit to the Arrhenius-type dependence N^i/N = Aexp(—Q/(K)) + B plotted with 
A = 0.12, B = 0.037, Q = 3.3 eV. The initial exponential drop in the number of dislocations is 
clearly seen at (if) -1 < 2 el/ -1 . Next, at higher (K) , the exponential decay is replaced by a 



power law decay 
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32] : in the figure N^isi /N ~ const due to the narrow range of the inverse mean 



kinetic energies involved. The simulation parameters were chosen to match the recrystallization 



experiment 
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First, the order parameter of the paired defects - dislocations, - was systematically lower 
for 7-fold cells. This is not surprising actually from a purely geometric point of view because 
the 5- fold cell in a pair is more compact. 

Second, simulations manifested that not only isolated pairs - dislocations (s 7 ), but also 
compact triplets like (5 7 5), quadruplets (5 7 7 5 ) etc., or even elongated defect chains were quite 
frequent. Actually they dominantly defined the symmetry of the entire particle suspension, 
't would certainly be promising to connect the cluster formation in ordered complex plasmas 
35 1 with the general percolation process known in many similar applications (see, e.g., 
36 
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Third, in such melted clusters, in agreement with recent experimental observations 
32], the defect density permanently decreased upon cooling. At higher temperatures in the 
beginning of the recrystallization process, while the mutual interparticle collisions were still 
frequent, the defect density dropped exponentially. Then, at lower temperatures, the decay 
rate significantly slowed down (see Figj2]). 

A sharp drop in the defect numbers followed by a quasi-saturation resembles the well- 



known situation [38[ in which both thermal activation and tunneling events occur. Hence, 
by analogy, the fact that in our case the system of defects behaves in a similar way could 
be naturally explained by an annihilation scenario which is presumably of the dissipative 



tunneling type 
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39| at lower mean kinetic energies. 



IV. COMPLEXITY OF NUCLEATION KINETICS 

Nucleation of dislocations is another important example of spontaneous symmetry break- 
ing on a scale of elementary cells. 

Whatever the melting scenario would be true, still there would remain a question what 
mechanism explains nucleation of the primary dislocation clusters. Recently this issue has 
been studied experimentally: Spontaneous nucleation of the edge-dislocation pairs (followed 
by their dissociation) has been successfully observed at the kinetic level in the experiments 



with plasma crystals 22j. Since the Burgers vector of the entire lattice is kept constant (e.g. 
zero) spontaneously created dislocations must be paired forming defect quadruplets of the 
type (5 7 7 5 ). (The Burgers vector characterizes the magnitude and direction of the crystalline 



lattice distortion by a dislocation 



40|.) These dislocation clusters were created in the lattice 
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locations where the internal shear stress exceeded a threshold. It has also been shown that 
even an elementary act of nucleation is in fact a multi-scale process consisting of the latent 
'pre-phase', prompt nucleation of a defect cluster, and dissociation of the cluster followed 



by the escape of free dislocations 
In the experiments 
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23|. 



23| it was suggested that the stress that finally caused nucleation 



was affected by the differential crystal rotation. The exact reason of nucleation, however, 
was difficult to determine consistently. In simulations the nucleation conditions are certainly 
easier to identify. 

To demonstrate nucleation in simulations a 'deformable' hexagonal cell is used (FigJS]). 
It confines a 2D cloud of equally charged particles interacting pairwise via the Yukawa (the 
screened Coulomb) force: 



1 %q + 



R 2 R 



(9) 



where = — Tj is the relative coordinate and R = |Rjj| is the distance between the 
particles i, j with the coordinates r^, r^; q is the particle charge and A is the screening length. 
The cell design is similar to that applied in 16j to simulate melting and recrystallization 



FIG. 3: Dislocation nucleation affected by a simple shear strain induced externally. The hexagonal 
simulation cell contains 631 interacting particles. The forces, deforming the cluster, are applied 
from left to right on top of the cell, and from right to left on bottom. The central horizontal 
part of the cell is kept non-deformed. The deformation rate is re = 0.05, where r is the Coulomb 



time-scale 



la], £ is the shear strain 



40l | . Three panels show the triangulated particle positions for 
the following time instants (from left to right): shortly before nucleation starts, immediately after 
nucleation, shortly after dissociation of the dislocation pair. The newly nucleated dislocations glide 



along the lines with the maximal stress. 
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process of the plasma crystal. 

The hexagonal simulation cell has the evident advantage of flexible shape, compared 
to, e.g., a parabolic cell confinement. Deforming the boundary of the cell, it is simple to 
manipulate the particles in a tractable way. An additional option of variable geometry 
enables an opportunity to separate or consolidate pure shear and simple shear deformation 



40J if desirable. The strain rate is controllable during deformation as well. 

Fig. [3] shows a simple-sheared particle lattice layer. At a properly chosen loading rate 
deformation affects the shear instability that ends up with nucleation of defect clusters 
in the bulk of the lattice layer. After a while, when deformation becomes stronger, the 
components of the clusters decoupled and the newly born free dislocations glided away in a 
similar manner as the dislocations observed in experiments. 

V. TOPOLOGY OF THE DISLOCATION CLUSTER 

Symmetry alternation is of primary importance for understanding nucleation of disloca- 
tion clusters. The compact cluster design is magic in the sense that the hexagonal symmetry 
of the particle system neatly turns into a nearly tetratic symmetry of the cluster core (like 
lead turns into gold when touched by the Philosophers Stone), see Fig. HI 

Despite an apparent simplicity of the cluster interior - only four nearest neighbor particles 
(marked ABCD in Fig. HJ), the centers of the 5- and 7-fold cells, are in the core, - to discover 
the cluster topology was certainly a challenge 41]. In our case the interparticle interaction 
potential is of the screened Yukawa type, hence more compact in contrast to the cx r~ 3 
interaction potential in case of magnetically interacting super-paramagnetic colloid particles 



41] . Thus there is a unique opportunity to verify whether the core topology 



considered in 
is universal. 

Let us start with a simple model treating the cluster as constituted of two point-like 
dislocations, which are set apart at a distance r and allowed to glide only along two fixed 
crystallographic planes separated by one lattice period a, so that r = a/ simp, where tp is 
the angle of mutual orientation of the cluster components with respect to the gliding plane. 
The interaction energy of the point dislocations having the counter-directed Burgers vectors 



is 



4lM44j: 



u M {r) = -^—Mc 2 tr [\n(-) + ^] + const. (10) 

7TV3 O T 
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FIG. 4: Topology of a four-component defect cluster, (black dots) Position of the centers of the 
5- and the 7-fold cells obtained in simulations. The simulation time is At/r = 120. The mean 
kinetic energy of the particles is kept constant at the level < K >= 7.3 ± 0.2 eV. The mean size 
of the cell edges is < s >= 1/4 < p >= 0.96 ± 0.04, where < p > is the mean cell perimeter. (As 
a scaling parameter the lattice constant is used.) Centers of the 5-fold (A, C) and the 7-fold (B, 
D) cells compose a nearly regular compact tetragonal structure. The gray quadrilaterals represent 



a few superimposed experimental clusters 
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23J shown for comparison, (inserts) Patterns of 



orientational preference revealed a visible 'quasi-hexagonal' trend in orientations for (left insert) 
the 5-fold and (right insert) the 7-fold components of the clusters. The mean radii are -R5 = 0.721 > 
Rj = 0.629. The bold and the dashed arrows indicate schematically the Burgers vectors of the 
individual dislocations. The components of the clusters (a) and (b) are shown in the bottom panels. 
The individual dislocations should have anti-parallel Burgers vectors since the Burgers vector of 
the entire defect group must be zero. This evidently results in the mirror-isomeric configurations 
of the escaping dislocations (the gliding planes are indicated by the dashed lines) . 
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It has a minimum, a stable ground state, at r/a — \/2 = 1/ simp. It corresponds to if = 45°, 
hence tetragonal symmetry of the cluster core might be considered as preferred. 

This prediction agrees noticeably well with the results of simulations of finite clusters: 
On average in Fig. [4] the edge-to-diagonal angle in the cluster core is < if >= 42° ± 2°. 

It is worth noting that the cluster core is nearly cyclic. A measure of it immediately 
follows from the famous Ptolemy's inequality valid for any quadrilateral: 

Pt = 8 i«» + > 1, (11) 

where si 4 denote the (ordered) sides, and are the diagonals of the quadrilateral. Over 
80 % of the recognized clusters have 1 < Pt < 1.03 for the simulation results shown in 
Fig. HI For comparison a hexagonal four-side cell corresponds to Pt hex = 2/y/3 = 1.1547. 

Note also that a stable defect cluster could not be obtained only by shifting positions of 
four central particles from a hexagonal configuration to tetragonal one. Such deformation 
would be reversible, hence unstable. A weakly deformed environment, impeding relaxation 
of the core particles back to the stable hexagonal configuration, is indeed a necessary 'lock' 
making the deformation plastic, i.e. irreversible. 

At sufficiently strong external stress even a stable cluster dissociates. Whatever is the 
orientation of the cluster as a whole, escaping dislocations can glide only along two crystal- 
lographic directions (along Burgers vectors, see Fig. H] (a, b)). This naturally explains the 



ity of the defect configurations revealed by the 



22| 



asymmetry of the escape directions and chira^ 
newly nucleated dislocations in experiments 

VI. SUMMARY 



Spontaneous symmetry breaking is a common and inherent feature of many systems in 
physics as well as other fields, it plays an important role, for example, from classical one- 
component plasmas to modern string representations 45j, from the evolution of the early 



universe [46] to the dynamics of a wide variety of small-scale systems |47| . Therefore it is not 
surprising that SSB is present also in the physics of plasma crystals. As example we have 
considered dislocations in plasma crystals which exhibit a spontaneous disordering, involved 
in the process of melting, and form clusters that are caused by a shear instability and show 
an interesting topological symmetry. 
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